Dynamics of Fluid Vesicles in Oscillatory Shear Flow 
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The dynamics of fluid vesicles in oscillatory shear flow was studied using differential equations of 
two variables: the Taylor deformation parameter and inclination angle 9. In a steady shear flow with 
a low viscosity rji n of internal fluid, the vesicles exhibit steady tank-treading motion with a constant 
inclination angle 9q. In the oscillatory flow with a low shear frequency, 9 oscillates between ±#o 
or around 9o for zero or finite mean shear rate j m , respectively. As shear frequency / 7 increases, 
the vesicle oscillation becomes delayed with respect to the shear oscillation, and the oscillation 
amplitude decreases. At high / 7 with j m = 0, another limit-cycle oscillation between 9o — n and 
—9o is found to appear. In the steady flow, 9 periodically rotates (tumbling) at high T]i a , and 9 
and the vesicle shape oscillate (swinging) at middle rji n and high shear rate. In the oscillatory flow, 
the coexistence of two or more limit-cycle oscillations can occur for low / 7 in these phases. For 
the vesicle with a fixed shape, the angle 9 rotates back to the original position after an oscillation 
period. However, it is found that a preferred angle can be induced by small thermal fluctuations. 

PACS numbers: 



I. INTRODUCTION 
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and synthetic cap- 
in flows has received 



growing attention experimentally, theoretically, and nu- 
merically in recent years. Under typical experimental 
conditions, these objects have a constant internal vol- 
ume V unless large shear breaks them. Lipid vesicles and 
RBCs essentially have a constant surface area A, because 
the number of lipids of their membranes is maintained 
constant in typical experimental time scales. Under these 
constraints, they show rich behaviors under flows. Flow 
can induce large shape deformations of these objects and 
lead to transitions of dynamic modes. 

An isolated fluid vesicle exhibits three types of dy- 
namic modes in a steady shear flow, with flow velocity 
v = ^ye x , where e x is a unit vector in the flow direction. 
When the viscosity of internal fluid r]- m and membrane 
viscosity ?y m b are low, the vesicle exhibit a tank-treading 
(TT) rotation with a stationary shape and a constant 
inclination angle 9 > [see Fig. Ufa)]. At high r]- m and 
7y m b, the vesicle exhibit a tumbling (TB) motion, where 
9 rotates [see Fig. Tib)}. Around the TT-TB transition 
viscosity with hi gh s hear rate 7, a swinging (SW) motion 
appears @, U & M EI EE EE HH HI, where 9 and the 
vesicle shape oscillate [see Fi g. W[ c)] . This motion is also 
called trembling @, H, IE EE HI] or vacillating- breathing 

, E3 • These three types of motion can be understood 
by the perturbation theories for quasi-spherical vesicles 
Gj, El, El] or a generalized Keller-Skalak (KS) theory 
for deformable ellipsoidal vesicles [2l[ . 
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RBCs and synthetic capsules also show the TT-TB 
transition with increasing viscosity rji n . The shear elastic- 
ity of their membranes gives an energy barrier to the TT 
membrane rotation (3o| . Since a sufficiently large shear is 
necessary to induce membrane rotation, transition from 
the TB to TT motions occurs with increasing 7 at low 
Vin f27l . l3o| . The synchronization and intermittency of 
9 and membrane rotations app ear at the transition re- 
gion between TB and TT [IE SI]- At much higher 7, 
this energy barrier becomes negligible, so the dynamics 
becomes similar to that of lipid vesicles. Thus, lipid vesi- 
cles, which have no shear elasticity, can be considered as 
a simple model system of RBCs in high shear rates. 

The relaxation dynamics of soft objects can be exam- 
ined using time-dependent flows. However, compared to 
steady-flow conditions, the dynamics in time-dependent 
flows has been much less investigated. Recently, a 
membrane wrinkling after inversion of an elongational 
flow (Io| and shape or orientational oscillation in struc- 
tured channels [5Cj were discovered for fluid vesicles. For 
RBCs, a shape oscillation under oscillatory shear flow 
was observed experimentally [32j and explained by the 
extended KS theory [4l|, which gives differential equa- 
tions of three variables, a shape parameter, the inclina- 
tion angle 9, and phase angle <f>. At high shear frequency, 
multiple limit cycles coexist. At middle shear amplitudes, 
RBCs show complex behavior such as intermittent oscil- 
lation and complicated domain boundary [4fj . 

In this paper, we study the dynamics of a fluid vesi- 
cle in oscillatory shear flow with 7 = j m + 70 sin(27r/ 7 t) 
using the generalized KS theory [2 if , which has two vari- 
ables, a shape parameter and the inclination angle 9. The 
reduction of variables from three to two allows us to in- 
vestigate the details of the dynamics. We investigate i) 
what kinds of oscillatory dynamics appear with respect 
to the TT, TB, and SW motions for zero mean shear 
rate (7 m = 0), ii) how the vesicle responds to the shear 
oscillation with finite mean shear rate 7 m , and iii) how 
the thermal fluctuations affect the dynamics. 
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The vesicle dynamics is described by several dimen- 
sionless quantities. The relative ratio of the volume 
V and surface area A is characterized by the reduced 
volume V* = W/{^Af/ 2 = {R v /R s ) 3 or the excess 
area A s = A/R% - 4tt = 4n{(l/V*) 2 / 3 - 1}, where 
Rv = (3F/47T) 1 / 3 and R s = {A/ An) 1 ' 2 . The relative vis- 
cosity of the inside fluid and membrane are = Tftn/^o 
and 7y* lb = r/ m b/r]oRs, where 770 is the viscosity of the out- 
side fluid. The shape relaxation time of the vesicles with 
bending rigidity k is given by r = 770-Rg/ft (for rj* n = 1). 
This time is used to define a reduced shear rate 7* = jr. 

The theory and results for the steady flow are ex- 
plained in Sees. Hi] and IIII( respectively. The results 
of the generalized KS theory arc compared with those 
of perturbation theories, experiments, and simulations. 
The dynamics for zero and finite mean shear rate is de- 
scribed in Sec. |W] and in Sec. E] respectively. In Sec. 
VH the effects of the thermal fluctuations are investi- 
gated. Summary is given in Sec. I VIII 



II. GENERALIZED KELLER-SKALAK THEORY 

Keller and Skalak [29[ analytically derived the equation 
of motion of vesicles or clastic capsules based on Jeffcry's 
theory [5l[. In the KS theory, the vesicles are assumed 
to have a fixed ellipsoidal shape, 



a-2 



(1) 



where denote the semi-axes of the ellipsoid, and the 
coordinate axes Xi point along its principal directions. 
The x\ and x 2 axes, with a\ > a 2 , are on the vorticity 
(xy) plane, and the x$ axis is in the vorticity (z) direc- 
tion. The Taylor deformation parameter is defined as 
Q'd = (Li - L 2 )/{Li + L 2 ) = (ai - a 2 )/(ai + a 2 ), where 
the maximum lengths in three directions are L\ = 2ai, 
L2 = 2a 2 , and L3 = 203. The velocity field on the mem- 
brane is assumed to be 



ai a-2 
— x 2 , — xi,0 
a,2 ai 



(2) 



The vesicle motion is derived from the energy balance 
between supply from the external fluid of the vesicle and 
dissipation inside the vesicle and on the membrane. 



dt 2 l 
B = /o(/i 



Bcos(26>)} 



/r 1 



1 + /afo£, - 1) + hhr,* mh 



(3) 
(4) 



The membrane- viscosity term was derived by Tran-Son- 
Tay et al. [30( . Factor B is a function of the vesicle shape 
and viscosity rj* n and rf^ h . Here, / , fx, f%, and f 3 are 
given by 

/o = 2/(01/03 + 02/ai) = i 1 - "d)/(1 + "d), 
fx = 0.5(ai/a 2 - a 2 /ai) = 2a D /(l - a%), 



(a) TT 




e J 








(c)SW ^ — 




I' I ( i . 1: (Color online) Sliced snapshots of fluid vesicles in 
steady shear flow obtained by numerical simulations at 
V* = 0.78 and 77* n = 1. (a) Tank-treading (TT) motion at 
? 7mb = an d 7* = 1-84 (n/ksT — 20), where the inclination 
angle 9 has a constant value, (b) Tumbling (TB) motion at 
r]* mh = 2.9 and 7* = 0.92 (/t/fasT = 40), where 6 rotates, (c) 
Swinging (SW) motion at r)^ b = 2.9 and 7* = 3.68 (k/Ic-bT = 
10), where 9 oscillates. The other parameters are the same as 
in Fig. 1 of Ref. [p. 



h = 
h = 

9 = 



0.5g(a 1 + ev 2 ), 
0.5E s R /(f 2 V), 



(a 2 + s)^ 2 (a 2 + s)-^ 2 (a 2 + sy^ds, 
ai/(aia 2 a 3 ) 1/3 , 



where E R is an integral of shear stress over the membrane 
surface j2fj|.[30j. 

The KS theory is extended to include the shape defor- 
mation of fluid vesicles on the basis of the perturbation 
theory [H M OS M of quasi-spherical vesicles in Ref. 
[2H ]. The equation of the shape parameter oly> is given 
by 



~dt 




A OF 

j*nV* dan 



A 1 sin(2<9) 
(5) 

where A = 45/8tt(32 + 23r/f n + 16rj^ h )V* and A x = 
30/ (32+2377^+1677^). Here, the terms of 77*^ are added 
in Aq and A\ based on the theory in Refs. jla, l!^. This 
revision improves the 77^ dependence of fluid vesicles in 
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Ref. [2l[. The first and second terms in the last paren- 
theses represent the forces of the bending elasticity of 
the membrane and the external shear stress, respectively. 
The free energy is given by F = (rc/2) J(d + C 2 ) 2 dA, 
where C\ and C 2 are the principal curvatures at each 
point of the membrane. Here F is numerically calcu- 
lated for ellipsoidal vesicles with (xi/ai) 2 + ix-ija^) 1 + 
(x^/a^) 2 = 1. The prolate (a\ > a 2 — 03) and oblate 
(ai = a 2 > 0.3) shapes are energy minima and maxima, 
respectively, and dF/dan diverges in the limit of maxi- 
mum extension, an — * a™ x (V*). Equations © and ([5]) 
are numerically integrated using the fourth-order Rungc- 
Kutta method with a time step At < 0.0005. 

For the data in Ref. [2lJ and in Sec. IIII1 the interpo- 
lations of F and prefactors in Eq.(|3]) are employed. For 
the dynamics in the oscillatory flow at V* = 0.9 and 
^mb = 0> the- fit functions are employed instead to avoid 
the artifacts of non-smooth functions: (l/n)dF/dar> = 
30a D - 820aj^ - 1.7/ N /ag ax - a D + 1.7/ yfiff*, h = 
0.3185+0.6af ) +14af ) +exp{124(a D -0.39)}, and ag ax = 
0.356. The difference in the phase boundaries calculated 
with the interpolation and fit functions is very small, 
A v * n < 0.04. 



III. VESICLE DYNAMICS IN STEADY SHEAR 
FLOW 

First, we briefly describe the vesicle dynamics in steady 
shear flow. Figure [2a ) shows the phase diagram with di- 
mensionless parameters S = (77t/3v / 3)7?7o-Rv 3 /kAs and 
A = y/Z A s /10tt(32 + 23?7* n + 16ry* lb )/24 calculated by the 
generalized KS theory; see also lines in Fig. E2a). The 
parameters S and A are introduced in the perturbation 
theory of quasi-spherical vesicles [l5| ■ 

At low shear rate 7* < 1 (S < 1), the vesicle keeps a 
prolate shape, and the original KS theory [Eq.([3])] gives 
very good predictions. For B > 1, a stable fixed point 
9 = 0.5 arccos(l/-B) exists, i.e., the TT motion occurs. 
For B < 1, there is no fixed point, and the angle 9 pe- 
riodically rotates (TB). As rj* or 77^ (A) increases, the 
transition from the TT to TB motion occurs, where B 
decreases from B > 1 to B < 1. 

On the other hand, higher shear at 7* > 1 (S > 1) 
can induce large shape deformation of the vesicles, and 
the SW phase appears between the TT and TB phases. 
In Eq. ([5]), the shear force depends on 9 as sin (29), so 
the shear increases ao (elongation) for < 9 < n/2, and 
decreases an (shrinkage) for — n/2 < 9 < 0. The SW 
motion is generated by this shape deformation as follows 
[2ll ]: A prolate vesicle starts 9 rotation with B < 1 like 
in the TB motion, and then the vesicle shrinks to more 
spherical shape at 9 < 0, which has greater B. When B 
becomes greater than l/cos(20), the right hand side of 
Eq.([3]) changes its sign and then the angle 9 increases. At 
9 > 0, the vesicle elongates to the prolate shape. Thus, 
9 and qd oscillate. 

The diagram of these three phases is also obtained in 
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FIG. 2: (Color online) Dynamics of fluid vesicles in steady 
shear flow, (a) Dynamical phase diagrams as a function of 
two dimensionless parameters S and A at V* — 0.85, 0.9, and 
0.95. The dashed lines represent the boundary of the TT 
phase. The solid lines represent the boundary between the 
TB and SW phases. In the upper three lines, rj* n is varied at 
^mb = 0. In the bottom three lines, 77^ is varied at r)* n = 1. 
(b) Stable angle 6 of the TT motion for rj* n — 1 and ?7^ b = 
as a function of V* (As in the inset). The solid (red) line rep- 
resent the KS theory [2!|. The dashed (blue), dotted (brown), 
and dashed-dotted (green) lines represent quasi-spherical per- 
turbation theories at 7 — > by Lebedev et al. [l5| and Danker 
et al. and at 7 — » 00 [Q, respectively. Th e sq uares and 
circles indicate the simulation results in Ref. [ill ] and Ref. 
[20I ]. respectively. The dotted (black) lines in the inset repre- 
sent a leading-order approximation of the perturbation theory 
by Seifert [IjJ (black), (c) Dependence of 9 and angular ve- 
locities on ?7* n at V* = 0.9 and rf^ h = obtained by the KS 
theory. 
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the perturbation theory of quasi-spherical vesicles [14l,[l5l 
[3 and by experiment Q. All of them give qualitatively 
the same diagram but quantitative differences exist. In 
the perturbation theory by Lebedev et al. [H, EH > the 
diagram is determined by only two parameters S and A. 
The TT-TB transition occurs at A = 2/V3 ~ 1.15 for 
5^0. The SW phase appears between the TT and 
TB phases at S > V3- The width of the SW phase 
is small, AA S w ~ 0.1 [1.4 < A sw < 1.5 at S = 10]. 
Danker et al. [14[ extended Lebedev's theory to take into 
account a higher order term. In Danker's theory, the SW 
phase becomes wider at larger As (1.4 < Agw < 1-7 at 
S = 10 and As = 1). Very recently, Deschamps et al. [1| 
obtained the phase diagram experimentally. Surprisingly, 
the resulting phase diagram depends only on A and S 
for the wide range of the excess area 0.2 < As < 2.2 
(0.98 > V* > 0.79) as Lebedev's theory However, the 
width of the SW phase is much wider, 1.5 < Asw < 2.2. 
In our generalized KS theory, the SW phase is shifted 
upwards with decreasing V* but the region of the SW 
phase is closer to the experimental results as compared 
to the perturbation theories; see Fig. [H^a). 

To compare the results in more detail, the stable an- 
gle 8 in the TT phase at r]* n = 1 and r) mh = is 
shown in Fig. Efb). At V* > 0.9, several theories 
[H 0Jll El, [221, simulations 0, Hfl H|, and experi- 
ments Ij, |6[ show good agreements. Kantsler and Stein- 
berg [6| reported that their experimental results agree 
very well with the leading order approximation of the the- 
ory ifi = tt/4 - C\/A"i) with C ~ 0.35 at V* > 0.82 
(As < 1.8); see the black dotted line in the inset of Fig. 
Mb). However, at the lower volumes V* < 0.8 (As > 2), 
deviations between the theories are seen. The KS the- 
ory predicts a gradual decrease in 8, while in the pertur- 
bation theories, 8 rapidly decreases and the TB motion 
occurs; see Fig. &b). The KS theory gives very good 
agreements with the simulation data of the boundary in- 
tegral simulation (without thermal fluctuations) [Til ] and 
multi-particle collision (MPC) dynamics (with thermal 
fluctuations) [2(| . Since the perturbation theory assumes 
As < 1, this TB phase at rj* n = 1 and 77. 



10F * * *~ 

(a) 



mb 



would 

be an artifact of the approximation. It is surprising that 
the perturbation gives good predictions even at As — 1. 
Thus, the perturbation theories should not be applied to 
the dynamics for As > 2. 

Previously, we simulated the SW motion in two- (Ref. 
J23|) and three-dimensional spaces (Ref. (2l|) at a re- 
duced area A* = 0.7 and volume V* = 0.78 using the 
MPC method [HH, respectively (see Fig. [1]). The region 
of the SW phase does not completely agree with the pre- 
diction of the generalized KS theory but it can be fitted 
by modifying the factors Aq or A\ . Since these factors are 
derived from the perturbation theory, it is less reliable at 
low reduced volumes. In this paper, we use the vesicle at 
V* = 0.9 and ?y* lb = to investigate the dynamics in the 
oscillatory flows, since this condition gives good agree- 
ments with experimental results. The rj* n dependence of 
8 and angular velocities in the KS theory (7* <C 1) arc 



1 r 



0.1 

0.15 
0.1 



0.05- 



X Alii * 



SW 1 



TB - 



— • • • 

4^ — I h 



!■ ■ ■ 

-H — + h 



- (b) 



y =0.1 




i a ! 1 



_l I L 



FIG. 3: (Color online) (a) Dynamical phase diagram at 
V* = 0.9, 77^ b = 0, and 7 m = 0. The closed circles, squares, 
and triangles indicate the TT, TB, and SW phases in the 
oscillatory flow. The dashed and solid lines represent the 
boundary of the TT, TB, and SW phases for the steady flow 
with 7* = 7o- (b) Bifurcation lines from one limit cycle to two 
limit cycles. The dashed-dotted (green), solid (red), dashed 
(blue) lines represent 70 = 0.1, 1, and 10, respectively. 



shown in Fig. [He). The mean squared angular velocity 
9m for the rotation is calculated as 8f n = (1/tt) f 8 d6. 
At the TT-TB transition point ?7* n = i]* = 5.3, deriva- 
tives of the rotational frequencies or mean angular ve- 
locities of the inclination angle 9 and phase angle u> arc 
discontinuous, while 9 m smoothly changes. 



IV. VESICLE DYNAMICS IN OSCILLATORY 
SHEAR FLOW WITH 7m = 

In this section, we describe the vesicle dynamics in 
the oscillatory shear flow with no net flow (j m = 0). 
Figures [3] and [4] show the phase diagram and domain 
boundaries of the phases. The dynamics in the oscillatory 
flow reflect the dynamics in the steady flow, but the phase 
boundary is smeared in the oscillatory flow. The SW and 
TB motions can occur at lower rj* n than that in the steady 
flow, since transient motions in the steady flow can be 
stabilized in the oscillatory flow. At low shear frequency 
/*, the vesicle repeats the motion for r)* n in the steady 
flow with 7* ~ 7q and 7* ~ — 7q. At a higher frequency 
/* , the vesicle shows two types of oscillations depending 
on the initial positions in the phase space. In the middle 
frequency /*, three or more limit-cycle oscillations can 
coexist at high r]* n (SW, TB). 

First, we show the dynamics at low r)* n . For a low 
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FIG. 4: (Color online) Domain boundary of the limit-cycle 
oscillations at (a) (70, T]* n ) = (1,1), (1,6), and (10,8); (b) 
(70, Via) = (b 4 ); and ( c ) (7o, Via) = (10,6). Each domain 
consists of the initial positions (cto,0) = (0.325, #i n t) at t = 
approaching the same attractor. 



shear frequency (/* < 0.1), the vesicle approaches one 
limit-cycle oscillation from any initial position, and the 
angle 9 stepwisely oscillates between ±0 O , where 9q is 
the stable angle of the TT phase in the steady flow; see 
Fig. [5jc). We call this oscillation as TT-based oscilla- 
tion. For a high shear frequency (/* > 0.1), the vesicle 
cannot relax to the angle ±0n of the steady flow for a 
half period l/2/ 7 . Instead, two limit cycles coexist as 
shown in Figs.[5£b) and (c): Oscillation between ±8q or 
between — 9q and #0 — 71". An approached limit cycle is 
chosen by an initial angle 9 but is almost independent 
of an initial shape ao- The trajectories on the («d, 9) 
plane in the stroboscopic map at t = n// 7 are shown in 
Fig. HJa), where n is an arbitrary integer. It rapidly re- 
laxes onto a sinusoidal curve and slowly approaches one 
of the stable fixed points. Two saddle points are seen 
on the curve at 9/tt = —0.18 and 0.38. Figure 2] shows 
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FIG. 5: (Color online) Vesicle dynamics in the oscillatory 
flow, (a)-(c) TT phase at (70, J7* n ) = (1, 1). (a) Stroboscopic 
map of the approach to stable fixed points [(blue and red) 
squares] from several initial positions on (ckd, 0) at t = n/f-y. 
The arrows represent the approaching directions. Time evo- 
lutions of QfD and 6 in the limit-cycle oscillations are given in 
(b) and (c), respectively, (d), (e) Trajectories of two limit cy- 
cles at high frequency /* = 0.25 for (70, f}ia) = (3, 1), (6, 10), 
and (8, 1). 



the domain boundary for initial angles 9- mt . Since the 
domain boundary is not sensitive to an initial value of 
an, the data with a constant initial shape = 0.325 
are shown. As /* decreases, the left stable and unstable 
fixed points become closer, and two domains merge at 
/* = 0.11; see Fig.QJa). Thus, the left fixed point (oscil- 
lation between —9q and #0 ~ it) at high /* is generated 
by a saddle-node bifurcation. The frequency /* for this 
bifurcation decreases with increasing r]* n ; see Fig. [H^b). 
This can be understood by slower relaxation velocity # m 
at higher r)* n ; see Fig. ^c). 

Figure [6] shows the dependence of the maxima and 
minima of cud and 9 in the TT-based oscillation. Here, 
we only show the data for < t < 0.5// 7 because 
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FIG. 6: (Color online) Dependence on shear frequency /* in 
the TT phase at rj* n = 1: (a) Maximum and minimum of od, 
(b) maximum of 8, and (c) time t at maximum and minimum 
of «d • The dashed (blue) and solid (red) lines represent 7* = 
1 and 10, respectively. 
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FIG. 7: (Color online) Time evolutions of qd and # in the 
limit-cycle oscillations for the TB phase at (70, T]* n )= (10,8). 
The solid (red) and dashed (blue) lines represent /* = 0.01 
and 0.021, respectively. 



FIG. 8: (Color online) Dependence on shear frequency /* 
in the (a), (b) SW and (c) TB phases at (70, r]*J)= (10,6) 
and (10,8), respectively, (a), (c) Maximum and minimum of 
c<D. The peak-to-peak amplitudes of 9 in the oscillations are 
shown in (b) and the inset of (c) . The solid (red) and dashed 
(blue) lines in (a) and (b) represent the data obtained with 
increasing or decreasing /*, respectively. 



a D (* + 0.5// 7 ) = a D (t) and 0(i + O.5// 7 ) = -6(t). In the 
low frequency limit /* < 1, the vesicles show maximum 
and minimum deformation at t = 0.25// 7 and t = 0, 
where the shear stress ?7o7(i) has maximum and min- 
imum values, respectively. As /* increases, the times 
at maximum and minimum deformation become delayed 
and approach t = 0.5// 7 and t = 0.25// 7 . On the other 
hand, always has maximum and minimum values at 
t = 0.5// 7 and t = 0, respectively. The amplitudes of qd 
and 9 in the oscillation rapidly decrease with increasing 
/* at /* > 0.1. This reduction of amplitudes and the 
coexistence of two limit cycles at high /* occur due to 
faster change of shear direction than vesicle relaxation. 

At low r]* n [(brown) circles in Fig. [21a)], one region of 
two limit cycles appears only at high /*. On the other 
hand, at higher r)* n [(blue) triangles and (green) squares 
in Fig. [3{a)] , regions of two or more limit cycles appear 
at lower /*. Next, we describe TB-based oscillation at 
a high T]* n [(green) squares in Fig. [3^a)]. At low /*, 
the vesicle rotates clockwisely at < t < 0.5// 7 and 
rotates back at 0.5// 7 < t < l// 7 ; see Fig. [71 Typi- 
cally, the shape of time-evolution curves are symmetric 
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FIG. 9: (Color online) Time evolutions of qd and 9 in the 
limit-cycle oscillations for the SW phase at (70, »7* n ) = (10,6) 
at (a)-(c) — 0.01 and (d) 0.02. Two coexisting limit cycles 
are shown as the solid (red) and dashed (blue) lines. 



like the TT-based oscillation, but asymmetric shape (see 
dashed lines in Fig. [7]) and the coexistence of two limit 
cycles can occur even at low /*. The rotation ampli- 
tude of 9 stepwisely decreases with increasing /* and the 
amplitude of «d oscillates as shown in Fig. Mc). At 
peaks of rotation amplitude in the inset of Fig. [8jc), two 
limit cycles with different rotation amplitudes are ob- 
tained with increasing and decreasing /*. Thus, narrow 
regions of two limit cycles repeatedly appear at /* < 0.1. 
At high /*, a wide region of two limit cycles appear like 
the TT-based oscillation. The trajectories of the oscil- 
lation around 0/tt = —0.5 are also very similar; see Fig. 
[SJe). However, the trajectories around 8/n = are con- 
vex upward and downward for the TB- and TT-based 
oscillations, respectively, (see Fig. GJd)), since 9 rotates 
in the opposite direction. 

At middle rj* n and high 7* [(blue) triangles in Fig. 
O^a)], the SW-based oscillation appears. The different 
frequencies of oscillations of cvd and 9 occur in a wider 
region of /* than the TB region; see Fig. The do- 
mains of multiple limit cycles have a complicated shape, 
as shown in Fig. life). These domains appear or disap- 
pear via saddle-node and Hopf bifurcations. The ampli- 
tudes of «d and 9 depend on trapped limit cycles. As /* 
increases or decreases, the vesicles are often trapped in 
different limit cycles; see Figs. IDJa) and (b). The peaks 
in Fig. [5{b) represent TB rotation, which coexists with 
SW oscillation. 




e n /n 



FIG. 10: (Color online) Vesicle dynamics for the TT phase 
in the oscillatory flow with finite mean shear rates -y m . (a-c) 
TT phase at r?* n = 1, 7^ = 2, and 70 /7m = 0.5. (a), (b) Time 
evolutions of an and 9 in the limit-cycle oscillations. The solid 
(red) and dashed (blue) lines represent — 0.01 and 0.2, 
respectively, (c) Trajectories at / 7 /7 m = 0.01, 0.1, 0.2, and 
0.5. (d) Return map sampled stroboscopically at t = n// 7 
in the TB attractor at r]* n = 8, 7 m = 10, and 70 /7m = 5. 
Circles (red) and squares (blue) represent ,f 7 /7 m = 0.1 and 
0.15, respectively. Solid (black) lines show 9 n+ i = 9 n — n and 
9 n +i — 9 n — 2ir. 



V. VESICLE DYNAMICS IN OSCILLATORY 
SHEAR FLOW WITH FINITE 7m 

Next, we consider the vesicle dynamics in the oscilla- 
tory flow with < 70 < 7m- In the TT phase, the shear 
rate 7 is always positive and the vesicle keeps a stable 
TT angle with 9 > 0. Unlike the oscillatory flow with 
7 m = 0, only one limit cycle exists even for high / 7 . The 
shape «d and angle 9 oscillate around the steady values 
at j(t) = 7 m ; sec Fig. [TU] Their amplitudes increase 
with increasing 70 and decrease with increasing / 7 ; see 
Fig. 111! In the low frequency limit, the maximum and 
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FIG. 11: (Color online) (a)-(c) Dependence on the frequency 
fj for the TT phase in the oscillatory flow with 7 m = 1 or 
2: Maximum and minimum of (a) ao and (b) 9; and (c) 
time t at maximum and minimum of an and 9. The solid 
(red), dashed-dotted (green), and dashed (blue) lines repre- 
sent (7m, 7o /7m) = (2,0.05), (2,0.25), and (1,0.25), respec- 
tively, (d) Shear-amplitude 70 dependence of the maximum 
and minimum of 9 at 7 m = 2. The solid (red), dashed-dotted 
(green), and dashed (blue) lines represent /* = 0.01, 0.1, and 
0.5, respectively. 



minimum of cud (minimum and maximum of 9) appear 
at t ~ 0.25// 7 and t = 0.75/ / 7 , where the shear stress 
?7o7(t) has maximum and minimum values, respectively. 
The minimum and maximum times are delayed with in- 
creasing / 7 , and these delays are almost independent of 
70 and 7 m ; see Fig. [TTTc). In the high frequency limit, the 
maximum and minimum of «d are reached at t = 0.5// 7 
and t = l// 7 , respectively, where 7 = 7 m . The angle 9 
shows greater delays than au , since stable 9 is varied not 
directly by 7* but by the evolution of au . The phase de- 



lay between au and 9 approaches 7r/4 as / 7 /7 m increases; 
see Fig. [TUTc). These oscillations also occur for RBCs, 
where the TT motion generates additional oscillation of 
the shape and angles [4l|. A similar phase delay is ob- 
served in the experiments of RBCs for oscillatory shear 

flow mi. 

In the TB and SW phases, the shear oscillation is 
weakly coupled with the oscillations of qd and 9. In most 
of the parameter regions, the TB or SW oscillation is not 
synchronized with the shear oscillation. The synchro- 
nization only occurs in very narrow regions, where the 
tumbling or swinging frequency has integer ratios with 
the shear frequency: /tb = { n / m )f~< °rfsw = {n/m)f<y, 
where m and n are arbitrary integers [531 ■ The synchro- 
nization with /tb = f-y and /tb = 2/ 7 appears only 
at 0.1522 < / 7 / 7m < 0.1542 and at 0.0764 < / 7 / 7m < 
0.0768, respectively, for r)* n = 8, 7^ = 10, and % = 5. 
Other synchronizations have much narrower regions. Fig- 
urc llOf d) shows the return map of typical attracted orbits 
in the TB phases. At 7o/7m = 0.15, an intermittent dy- 
namics appears; 9 n+ i (6 at t = (n + l)// 7 ) approaches 
d n +i = n — 7r at 9 n = (but not contacted), so that the 
vesicle is trapped on 9 = at t = nj / 7 and intermittently 
escapes. 

Elastic capsules such as RBCs [27j and synthetic cap- 
sules [ID, [43l 0, HH, show similar oscillations of au 
and 9 at the tank-treading phase in steady shear flow. 
To overcome the energy barrier, the TT membrane ve- 
locity w of the elastic capsules oscillates with the TT ro- 
tation frequency, and induces the oscillations of qd and 
9 [H, |4Q]. For the elastic capsules, the phase difference 
between au and 9 is also 7r/4 at a large shear rate, but 
approaches at lower bound j tt of the TT phase [4C| . 
This phase shift occurs since 9 oscillation is caused by 
tu oscillation as well as q;d oscillation unlike for the lipid 
vesicles in the oscillatory flow. 



VI. EFFECTS OF THERMAL FLUCTUATIONS 
IN OSCILLATORY SHEAR FLOW 

In the small shear limit 7* — + 0, the vesicle deforma- 
tion becomes negligible and the vesicle dynamics is given 
by the original KS theory, i.e., a^> and the factor B in 
Eq. are constant. The original KS theory gives no 
attractor in the oscillatory flow with 7^ = 0. Since Eq. 

is symmetric to time reversal, 9 exactly returns to the 
original angle at t' = t + n/ / 7 ; see Fig. [TSta). The shape 
deformation given by Eq. ([5]) breaks this reversibility 
but the approach to the limit cycle is slower at smaller 
7p. At small Jq , the thermal noises are typically not 
negligible. Recently, Pine et al. [54[ demonstrated that 
noises can break the time reversibility: The interactions 
between spherical particles can induce the diffusion of 
particle positions at t = n/ f-y in the oscillatory shear 
flow. In this section, we clarify the thermal noise effects 
on the vesicle dynamics with a fixed shape. 

Gaussian white noise g(t) is added to Eq. ([3]) to take 
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into account the thermal fluctuations, where (g(t)) = 0, 
(g(t)g{t')) =2D5{t-t'). The fluctuation-dissipation the- 
orem gives the diffusion constant D = /cbT/C where ( 
is a rotational friction coefficient and k^T is the thermal 
energy; for a sphere, C, = 87r?7oi?g. A dimensionless quan- 
tity, the rotational Peclet number \ = 70/-D = ^Q/ksT 
represents the shear amplitude relative to the thermal 
fluctuations. Equation (J3]) is numerically integrated us- 
ing the second-order Runge-Kutta method [55[ with a 
time step At < 0.00005. 



First, we consider the TT vesicle at r)* n = 1 (B = 1.7, 
where the stable TT angle 9q = 0.1 57r in the steady flow). 
At x iS 1? the thermal fluctuations are dominant, and 9 
show a uniform distribution. As x increases, a peak grows 
at 9 ~ —6*o or at 9 ~ 9q in the distribution of at t = 
n/f 7 or at t = (n + 0.5)// 7 , respectively; see Figs. [T^b) 
and (c). For n// 7 < t < (n+0.5)// 7 , a half of the vesicles 
rotate 6 ~ — 9q to 9q and the other half rotate to 9q — it. 
Thus, the peak at 9 ~ — 9 splits to two peaks, and then 
they merge into a new peak at 9 ~ 9o- Thus, the thermal 
fluctuations lead the vesicle to the orbit of the deformed 
vesicle; compare Fig. EJc). These effects can be explained 
by the orbits in Fig. [T27 a). The angles in the wide region 
of —0.8 < 9/tt < 0.15 at t = (n + 0.5)// 7 accumulate in 
the narrow region —0.2 < 9/tt < 0.1 at t = (n + l)// 7 . 
The combination of this accumulation and the thermal 
diffusion gives the peaks in the angular distribution. The 
width of the distribution is calculated as the standard 
deviation xj {(A9 / V) 2 ) for 6» max -0.57r < 9 < 6> max + 0.57r, 
where (9 max is the maximum angle; see Fig. IT2T d). For 
the uniform distribution, v /((A9/tt) 2 ) = l/2y/3 ~ 0.29. 
At x ^ 1; t nc width decreases with decreasing /*. A 
sharper peak is obtained for smaller /* with large x- 

Next, we investigate the vesicle in the TB phase. In 
the limit r)l n — > 00 , the vesicle becomes a solid object, 
where no shape deformation occur even at high 7* . Then, 
the dynamics is given by Jeffery's theory [5l|, where 
B = /0/1 in Eq. ([3j> (B = 0.59 for a prolate object at 
V* = 0.9). In the absence of the thermal noise, the orbit 
is time reversible and the angular distribution has a peak 
at 9 = in the steady flow. The thermal fluctuations 
slightly shift this peak to positive 9 [20], and the peak 
becomes smaller at smaller x\ see Fig- IT5Tc). In the oscil- 
latory flow, the peak in the angular distribution moves to 
negative (positive) 9 for t = n/f^ (t = (n + 0.5)// 7 ) with 
increasing /*; see Fig. [T31 For n// 7 < t < (n + 0.5)// 7 , 
most of the TB vesicles show clockwise rotation from 
9 ~ — 7r/4 to — 37r/4, while the peak moves from — 7r/4 
to 7r/4. The sharpest peak is obtained for /* ~ 0.1 with 
large x m contract to a monotonic dependence on /* in 
TT phase; compare Figs. [T2Td) and flUT c). At smaller 
frequency /* <C 0.1, it approach the uniform distribu- 
tion. Thus, very small thermal noise can break the time 
reversibility for both the TT and TB phases and generate 
a preferred orbit, although the relaxation time is longer 
for the smaller noise. 
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FIG. 12: (Color online) Dynamics of vesicles with a fixed 
prolate shape at V* = 0.9 and r)* n = 1 in the oscillatory 
flow with -y m = 0. (a) Time evolution of 6 without thermal 
fluctuations for /* = 0.1. (b) Time evolution of probability 
distribution of at Peclet number \ = 100 an d fj — 0.1. (c) 
Probability distribution of at t = n// 7 for various \ with 
/* — 0.1. (d) Peclet number x dependence of the standard 
deviation of 9 distribution at t — n/f^ for various /*. The 
error bars are smaller than symbols. 



VII. SUMMARY 



We found the fruitful dynamics of fluid vesicles in os- 
cillatory shear flow. For no net flow (7 m = 0), the vesicle 
approaches one or multiple limit-cycle oscillations. In the 
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FIG. 13: (Color online) Dynamics of a solid prolate object 
(?7* n — > oo) at V = 0.9 in the oscillatory flow with 7 m = 0. 
(a) Time evolution of probability distribution of 6 at Peclet 
number \ = 1000 and /~ = 0.1. (b) Probability distribution 
of at t = n// 7l x = 1000, and /* = 0.05,0.1,0.2,0.5. (c) 
Peclet number \ dependence of the standard deviation of 6 
distribution at t — n/f 1 for various f^. The data (A) for 
steady flow with \ = j/D are also shown. The error bars are 
smaller than symbols. 

TT phase (low inertial viscosity r]* a ) at a low shear fre- 
quency /*, the inclination angle 9 oscillates between the 
steady TT angle ±9q in the steady flow with 7* ~ ±7q 
at n < fjt < n + 0.5 and n + 0.5 < / 7 i < n + 1, respec- 
tively. In the TB phase (high r]* n ), 9 rotates clockwisely 



and then rotates back. In the SW phase (middle r/* n and 
high 7q), the vesicle shape and 9 oscillate. At high /*, 
two limit cycles coexist. The angle 9 oscillates around 
9 = or 9 = — 7r/2. At some regions of middle /* two 
or more limit cycles can exist in the TB and SW phases. 

In the oscillatory flow with the finite net flow at 
7m > 7o > 0, only one limit cycle exists in the TT 
phase. As /* increases, the oscillation of shape and 9 
become delayed compared to the shear oscillation. The 
SW and TB oscillations arc typically not locked with the 
shear frequency. In the limit of small 7q or high rj* n (cor- 
responding to solid objects), the dynamics is integrable 
and no attractor exists in the absence of noise. How- 
ever, we found that the thermal fluctuations can induce 
an attractive orbit. 

In the previous paper [Z(| , we investigated the dynam- 
ics of an RBC in oscillatory shear flow using the differ- 
ential equations of three variables. RBCs show TT, TB, 
and intermediate motions depending on 7g and r]* n . In 
the TT and TB phases, the dynamics of RBCs is very 
similar to that of fluid vesicles, while the membrane en- 
ergy barrier of RBCs induces additional shape and an- 
gular oscillations. For RBCs with the intermediate shear 
amplitude, the small variations in the parameters can 
induce the dynamic modes from intermittent to synchro- 
nized oscillations, and change the number of coexisting 
limit-cycle oscillations. 

Recently, Deschamps et al. demonstrated that a 
general flow controlled by a 4-roll mill device is a good 
experimental setup to investigate the details of the vesi- 
cle dynamics. Our study shows that the oscillatory shear 
flow is another condition to reveal the dynamics of vesi- 
cles, cells, and capsules. The analysis of the frequency de- 
pendence helps understand their viscoelastic behaviors. 
We hope that our study stimulates further experimental 
study on oscillatory flows. 
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